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Abstract 

The Lorentz-Dirac radiation reaction formula predicts that the position shift 
of a charged particle due to the radiation reaction is of first order in acceler- 
ation if it undergoes a small acceleration. A semi-classical calculation shows 
that this is impossible at least if the acceleration is due to a time-independent 
potential. Thus, the Lorentz-Dirac formula gives an incorrect classical limit 
in this situation. The correct classical limit of the position shift at the lowest 
order in acceleration is obtained by assuming that the energy loss at each 
time is given by the Larmor formula. 
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A classical charge e with four-velocity emits electromagnetic radiation with energy- 
momentum P M according to the Larmor formula: 



where r is the proper time along the world line of the charge e and where a = ^J—u^u^ 
(with ii^ 1 = du^/dr) is the acceleration of the charge. (See, e.g., Ref. [|1].) We have defined 
a = e 2 /47r. It is widely believed that the radiation reaction is given by the Lorentz-Dirac 
formula The corresponding four- force is given by 

tf, = _(- V+-^J. (2) 

One can readily verify that u^K^ = 0. (This is necessary for consistency of Newton's 
equation mii^ = K^, where m is the mass of the charge.) Moreover, if the acceleration lasts 
only for a finite time, the total energy-momentum lost due to this four-force agrees with the 
total energy-momentum radiated in the form of electromagnetic radiation according to the 
Larmor formula (|l|). 

Although the Lorentz-Dirac formula (fj) is widely accepted, it is known to have some 
unphysical features. (See, e.g., Refs. |]|4[].) Let us consider one-dimensional motion and let 
u tj. = (^o^ 1 ) — (cosh/3, sinh/3). Then the Lorentz-Dirac force is = |a(sinh/3, cosh/3)/3. 
Newton's equation reads 

d(3 2ad 2 (3 

where F ext (r) represents an external force. This equation allows runaway solutions. For 
example, if F ext (r) = 0, then f3{r) = /3 exp(3mr/2a) is a solution. Runaway solutions can 
be excluded by rewriting Eq. (|) as follows: 

d/3 



m— = e~ s F ext (r + r s)ds, 
dr Jo 



where tq = 2a /3m. However, this equation violates causality. This causality violation 
is usually dismissed as unimportant [|I] because its time scale tq is very small. (It is of 
order 10~ 24 sec for electrons.) This nevertheless casts some doubts on the validity of the 
Lorentz-Dirac theory. In this Letter we will show that in fact the Lorentz-Dirac theory does 
not describe the classical limit of the quantum theory if the acceleration is due to a static 
potential varying in one dimension, say in the ^-direction. (Moniz and Sharp concluded that 
the Lorentz-Dirac formula is reproduced in the classical limit by studying the Heisenberg 
equations ||. It is difficult to find the reason why we disagree since our method is quite 
different from theirs.) 

Once we start considering a quantum system with a static potential, a natural alternative 
classical theory emerges. The z-component of the momentum is not conserved in the rest 
frame of the potential. However, the energy is conserved. Hence, it is natural to expect 
that the correct classical limit is obtained by assuming that the energy is lost at each time 
according to the Larmor formula ([!]), disregarding the momentum conservation in the in- 
direction. We call this theory the Larmor theory in this Letter. We show that the Larmor 
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theory is likely to be the correct classical theory. We use natural units h = c = 1 unless 
otherwise stated throughout this Letter. 

We consider a charged particle moving in the ^-direction under the influence of a potential 
V(z) which is constant if \z\ is sufficiently large. (Thus, the acceleration lasts only for a 
finite time.) Let us first clarify how the Lorentz-Dirac and Larmor theories differ. 

We work in the non-relativistic approximation. Let us assume that the kinetic energy is 
dominant so that the particle always moves in the positive z-direction. The energy is defined 
as E = |mt> 2 + V(z), where v is the velocity of the particle. The rate of energy change in 
the Lorentz-Dirac theory reads 

dE D 2 .. 

— - — = -aw , (4) 
dt 3 v ; 

where v = d 2 v/dt 2 . On the other hand, according to the Larmor theory this quantity is 

dE L 2 . 2 

To first order in a, one can let v on the right-hand sides of Eqs. (§) and ([5]) be the velocity 
without any radiation reaction. By imposing the condition that Ed = El for t — ► — oo one 
has 

Note here that this quantity vanishes for t — > +oo because V(z) is assumed to become 
constant for large \z\. Thus, the two theories give the same final velocity. 

Now, let z denote the position without the radiation reaction. Let Adv and Adz be the 
corrections to the velocity and position at the lowest order in a according to the Lorentz- 
Dirac theory and A^v and A^z be the corresponding quantities for the Larmor theory. 
Then, using V'(z) = —mi), A D v = (d/dt)(A D z) and A L v = (d/dt)(A L z), one finds from 
Eq. ® 

d / A D z — A L z\ 2a d 



m— = -lnt>. 

dt \ v J 3 dt 

Let the initial and final velocities without the radiation reaction be Vi and Vf, respectively. 
Then, for t — > +oo, one has 

6m Vi 

Now, let a(t) be the acceleration. Then, to lowest order in a(t) one finds 

2a 2a 
A D z-A L z = —(v f -Vi) = — a(t)dt. (7) 

oTTl oTfl J— oo 

Thus, the two theories differ in the position shift, and the difference is of first order in 
acceleration. 

Next let us derive the formula for the position shift at t = for the Larmor theory, 
assuming that a(t) = for t > for simplicity. The energy formula (^) gives 

d i . . . . 2a , , 9 
-(vA L v-vA L z) = --a(tf. 
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From this we find the position shift at t = to lowest order in a(t) as 

(A L z) = ~ f° ( f ait'fdA dt 

Sp J-oo \J-oo J 

= — / ta(t) 2 dt, 

3p J-oo 

where we have replaced mv by the average momentum p. Define d(k) to be the Fourier 
transform of a(t), i.e., 



Then we find 



n+OO 

a(k) = I dta(t)e iht . 



2ia r+°° dk ■ 



Note that this position shift is of second order in acceleration. Since the difference in the 
position shift for t — > +00 given by Eq. ([/]), which equals that at t = if a(t) = for t > 0, 
is of first order, the position shift (Ajjz)q at t = according to the Lorentz-Dirac theory 
is of first order in acceleration. It will be shown next that the position shift is of second 
order in acceleration in quantum theory and, therefore, that the Lorentz-Dirac theory gives 
an incorrect classical limit. 

The quantum system corresponding to the classical one considered above is given, to 
lowest non-trivial order in e, by the following Hamiltonian: 

* 2 

H = |- + V(z) + H cm - -A ■ p, (9) 
2m m 

where H em is the Hamiltonian describing the free electromagnetic field and where p = 
—id/dx.. The relevant Hilbert space is the tensor product of the space of non-relativistic 
one-particle Schrodinger wave functions and the Fock space of photons. The underlying 
theory is the scalar QED with an external potential for the charged scalar field. (See, e.g., 
Ref. for the non-relativistic version of this theory.) One arrives at the Hamiltonian (JD 
by adopting the Coulomb gauge and dropping the terms of second order in e.[] 

First let us consider wave functions evolved in time by the Hamiltonian H = p 2 /2m + 
V(z). They represent accelerated particles without radiation reaction. Let V(z) = V = 
const, for sufficiently large and negative values of z and V(z) =0 for sufficiently large and 
positive values of z. Writing the solution to the time-independent Schrodinger equation as 
<f>p(z) exp(iP x x + iPyy), where the energy is E = (P 2 + P 2 + P 2 )/2m, one finds 



1 d2 t /y \ 
2^d? + Viz) 



P 2 



l The term proportional to e 2 <fi contributes only to the mass renormalization at the lowest 
order. The term representing the Coulomb interaction produces a contribution corresponding to 
Coulomb scattering of the charged particle off a vacuum loop under the influence of V(z). This 
contribution vanishes if one assumes that the potential V{z) is the same (including the sign) for 
both the positively and negatively charged particles. We will make this assumption here. 
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Let us assume that P 2 /2m 3> \V(z)\ and that V(z) is smooth enough so that the WKB 
approximation is applicable. The right-moving solutions (j)p(z) are given by 



m*) = «(*r 1/2 



cxp 



k(z) dz 



where k{z) = [P 2 - 2mV(z)] 1 / 2 . 

Let the initial state be ^ — "00 <S> 1 0) , where |0) is the vacuum state of the electromagnetic 
field and ip is a wave packet of the charged particle given by 



dPdP 



; f(P,P_ L )VPMz)e 



D iP ± - X± -iE(P,P ± )t 



1. 



Here, P ± = (P x ,P y ), E(P,P ± ) = (P 2 + P 2 ± )/2m, x ± = (x,y) and / dPd 2 P ± \f(P, P ± )\ 
The function f(P, P±) is assumed to be non-zero only for the values of P much larger than 
\2mV(z)\ l l 2 . The wave packet is assumed to be sharply peaked about a classical trajectory 
near the region with V'(z) ^ 0. The electromagnetic field A(x) in the Schrodinger picture 
can be expanded as 



a(x) = £ 



d 3 k 



(2tt) 3 2A; 



ikx 



ikx 



where are the (real) polarizations of the photon satisfying • k = and where k — |k|. 
The annihilation and creation operators a k j and a^ - satisfy the usual commutation relations 

[a k 'i'> a kj] = fyj<K k ' - k ) and [ak'j')°kj] = l a k>j>i a kj} = 0- 

By straightforward application of the time-dependent perturbation theory with the in- 
teraction Hamiltonian Hi = — (e/m)A ■ p, one finds the first-order component in the final 
state as 



d 3 k 



-J(P,P ± )\pp ± kj) 



x J dz'Mz 1 )^ ■ P ± (t> P (z') + Mz')¥ kzZ \ 



where 



\ PP± kj) = ^/p-Mzy p± ' x± ~ tE(p ' p±)t 



4 ? |o), 



p ± + k ± and P = Jp 2 + p 2 L -P 2 L + 2mk. 



At the lowest order in the WKB approximation, one has 4>'p(z) ~ i[P 2 — 2mV(z)] 



V2, 



>P [Z 



By choosing the wave packet such that |Pj_| <C P if f(P,P±) ^ 0, the term proportional 
to ■ Pj_ can be neglected. One can then choose the two polarizations so that = 0. 
Thus, one has 



. r d 2 p±dp r 



d 6 k 



e^I pk f(P,P)\pp ± kj), 



where 
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Ipk = M / dz Jp 2 - 2mV(z) <f) p (z)Mz)e- ikzZ 

V r J — oo 

/+oo r rz 

exp i / g(z')dz' — ik z z . 
-oo L JO 

We have defined 

g (z) = [P 2 - 2mV(z)} 1 / 2 - [p 2 - 2mV(z)} 1/2 . (10) 

It is assumed here that the WKB approximation is applicable to the final-state wave func- 
tions as well. This assumption is valid unless the photon carries away most of the kinetic 
energy of the charged particle. The phase space for such final states can be made arbitrarily 
insignificant by increasing P. Hence, it is reasonable to assume that the contribution from 
such final states can be neglected. 

The integral I p k can be rewritten by integrating by parts and dropping the surface terms 

as 



' pk 



-°° g\z) 

dz -r-—: -t-tt; exp 



oo 







2 

g{z')dz' 



W) - K\ 2 

(The surface terms can be dropped because we are dealing with a wave packet here, and 
the final result would be the same if we introduced a damping factor which made the sur- 
face terms vanish.) Now, the function g{z) in Eq. ( [10|) is multiplied by ft -1 if one restores 
h. Hence, in the classical limit % — > 0, the integral I p k tends to zero due to rapid oscilla- 
tions if P p. Since we are interested in the classical limit, we make the approximation 
P »s p. Now, by using conservation of transverse momentum, Pj_ = p± + k^, one finds 
|Pj_ — pj_| < fc(|Pj_| + | pj_ | ) - Note that m 3> |P_J + |p±| because the charged particle 
is assumed to be non-relativistic. Therefore, we conclude that |Pj_ — pj_| 2mk. From 



P = yp 2 + Pj_ — Pj_ + 2mk we have P ~ [p 2 + 2mk] 1 ^ 2 . Thus, in the non-relativistic and 
classical regime, one can make the following approximation: 

g (z) = [P 2 - 2mV(z)} 1/2 - [p 2 - 2mV(z)} 1/2 



V;, 



where v p (z) = [p 2 — 2mV(z)] 1 ^ 2 /m is the velocity of the classical particle with final momen- 
tum p in the z-direction. By using this approximation one finds 



rOO 



kv' (z)e 



-ik z z 



'I dZ [k-k z v p (zW eXP 



h ' -dz' 



Jo v p (z') 

Noting that k ^> k z v p (z) because the particle is non-relativistic, one finds 

^ r+oo ^ 

I P k~ T dta p (t)e lkt = T a p (k), 
k J-oo k 

where a p (t) is the acceleration of the classical particle with the final momentum p, with 
t = at z = 0. Thus, the first-order part of the final-state wave function in the momentum 
representation is 

'(2ir) 3 2k 3 
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We have let P ~ p and Pj_ ~ pj_. The probability of emission is obtained by integrating 
\i/ji\ 2 over p, p± and k. Here, we note that the amplitude of emission is proportional to 
the acceleration. Hence, the position shift must be at least of second order in acceleration 
because the probability is proportional to the amplitude squared. Thus, the Lorentz-Dirac 
theory cannot give the correct classical limit. 

By recalling that / dpd 2 p±\f(p, p_i_)| 2 = 1, the emission probability is found to be 



_ 2e 2 f d 3 k o 4a f°° dk |A .... 
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where p is taken to be the average value for the wave packet. We have used the fact that 
the average of (e^) 2 is 2/3. The probability V is infrared divergent if a p (0) ^ 0, i.e., if the 
initial and final velocities are different. The expected energy of the emitted photon is 

e 2 f 00 dk.„ ,, , „ 2a f+°° , , , 9 

671 Jo llX 6 J-oo 

Thus, the classical Larmor formula (|I|) is reproduced. 

Now, the position operator (in the momentum space) of the final-state wave packet 
evolved back to t — can be approximated by id/ dp. (The time t = is automatically 
picked out because the phase factor describing the time dependence of the wave function is 
chosen to be one at t = 0.) Let us concentrate on the sector with fixed and suppress p_i_. 
The contribution to the expectation value of id /dp from the tree diagram of first order in e 
is given by 

(z) = J rf 3 k(^,^) 
4e 2 r d 3 k 



3 J (2Tr)32kZ { ^ k)f ( p), ^ k)f ( p)) - 



where 



(A,B)= l -Jdp 



^dB_ _ dA_ B 

dp dp 



The acceleration is given by —V'(z)/m as a function of z, and its p and t dependence is only 
through z =pt/m. Thus, one has [t(d/dt) — p(d / dp)]a p (t) = 0. Hence, 

^-a p (k) = --d p (k) - -^-d p (k) . 
dp p pdk 



By using this formula and the fact that d p (—k) = d p (k) (because a p (t) is real), one obtains 

An r°° dk 

(*)o = y</(p),/(p)> I !^M^)I 2 

2ia f+^dp^, 2 (+°° dk - n ^ d 



The first term is the expectation value of the operator id /dp (without radiation reaction) 
times the probability of emission. This term will be cancelled by the "one-loop" correction 
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without emission at this order due to unitarity. This means that the expectation value of 
the final-state position shift evolved back to t = is given by 



/A v 2ia [ +00 dp. £ . 2 r+°°dk——d A 

It is clear that one obtains Eq. (pi) in the limit where the wave packet becomes concentrated 
at a point in the momentum space. Thus, the Larmor theory gives the correct position shift. 

In summary, it was shown in this Letter that the Lorentz-Dirac radiation reaction formula 
does not give the correct classical limit of the quantum mechanical system with a static 
potential, by demonstrating that the position shift due to that formula is incorrect. It was 
also shown that the correct classical limit of the position shift is obtained by assuming that 
the energy is lost according to the Larmor formula, disregarding momentum conservation. 
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